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Abstract 

X 

^ ■ We study the production of a large-p^ photon in association with a jet 



in proton-proton colhsions. We examine the sensitivity of the jet rapidity 
distribution to the gluon distribution function in the proton. We then assess 
the sensitivity of various photon -|- jet correlation observables to the photon 
fragmentation functions. We argue that RHIC data on photon-jet correlations 
can be used to constrain the photon fragmentation functions in a region which 
was barely accessible in LEP experiments. 
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1 Introduction 



The phenomenology of prompt photons is very rich and interesting, as the photon 
on one hand can be considered as a pointhke particle described by QED, leading to 
clean experimental signatures. On the other hand, the photon is also involved in 
hadronic phenomena, like the fragmentation of an energetic parton into a large-pr 
photon and hadronic energy. 

At the LHC, photon+jet final states will be important for jet calibration and pdf 
studies [H [21 [3l m O [6] . Diphotons will play an important role in the search for a 
Higgs boson with mass below ~ 140 GeV, where the decay into two photons is a very 
prominent channel [7l[8], for which the branching ratio is small, ~ O{10~^), and the 
signature is provided by a narrow peak over a huge background involving various 
components. Besides the so-called irreducible background from prompt diphotons, 
the background called reducible comes from photon-jet and jet-jet events, with the 
jet faking a photon in various possible ways (high px vr^ or other neutral hadrons, 
charged particles inducing the radiation of energetic photons e.g. by bremsstrahlung 
due to interactions with innermost layers of the detector, etc.). An accurate knowl- 
edge of the photon+jet rate in particular is required to estimate and control the 
reducible background to the Higgs boson search in the diphoton channel [9]. In 
addition, highly energetic photons are important signatures for various scenarios of 
physics beyond the Standard Model. Therefore, issues like controlling photon iso- 
lation or further constraining the parton-to-photon fragmentation functions are of 
major importance, and data from RHIC and the Tevatron should be exploited as 
much as possible to this aim. 

The production of prompt photons in hadronic collisions may be schematically 
seen as originating from either of two mechanisms. In the first one, which may be 
called 'direct' (D), the photon behaves as a high px colourless parton, i.e. it takes 
part in the hard subprocess, and is most likely well separated from any hadronic 
environment. In the other one, which may be called 'fragmentation' (F), the photon 
behaves hadron-like, i.e. it results from the collinear fragmentation of a coloured 
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high pt parton. In the latter case, it is most probably accompanied by hadrons - 
unless the photon carries away most of the transverse momentum of the fragmenting 
parton. 

From a technical point of view, (F) emerges from the calculation of the higher 
order corrections in the perturbative expansion in powers of the strong coupling a^. 
At higher orders, final state coUinear singularities appear in any subprocess where a 
high Pt outgoing parton of species k (quark or gluon) undergoes a cascade of succes- 
sive coUinear sphttings together with the coUinear emission of a photon. The higher 
order corrections to the cross section can be split into 1) a contribution free from 
these final state coUinear singularities, to be added to the Born term so as to build 
(D), and 2) a contribution (F) involving these singularities together with accompa- 
nying large coUinear logarithms. In (F), the final state coUinear singularities and 
accompanying logarithms can be factorised to all orders in from short distance 
terms according to the factorisation theorem, and absorbed into fragmentation func- 
tions of parton A; to a photon D^/k{z, M^) . Let us mention however that the splitting 
of the cross section between (D) and (F) is not unique and that the Djfk{z, M'^) 
depend on the arbitrary factorization scheme specifying which non-singular parts 
are factorised together with the coUinear singularities; the latter depend in partic- 
ular on some arbitrary fragmentation scale M^. We therefore need to define the 
scheme used. In this article, (D) is defined as the Born term plus the fraction of the 
higher order corrections from which final state coUinear singularities and accompa- 
nying coUinear logarithms have been subtracted according to the MS factorisation 
scheme. (F) is defined as the contribution involving a fragmentation function of any 
parton into a photon defined in the MS scheme. The partonic cross section can thus 
be written schematically as: 

da^ = d(7(^)(/x^M^M|)+ Yl dafH/i^M^M|)®D^/fc(M|) (1.1) 

k=q,q,g 

where /i, M, Mp are respectively the (arbitrary) renormalization, initial state fac- 
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torisation and final state fragmentation scales, and "®" stands for a convolution 
over the fragmentation variable. The point-like coupling of the photon to quarks 
is responsible for the well-known anomalous behaviour of D^/fc(z,M^), roughly as 
aem/ otsi^l)-, "^hen the fragmentation scale M^, chosen of the order of a hard scale 
of the subprocess, is large compared to (9(1 GeV). More generally, while the Mp evo- 
lution of these fragmentation functions is given by inhomogeneous evolution equa- 
tions whose kernels are computable in perturbative QCD, the z profiles of these 
fragmentation functions are not fully predictable from perturbative QCD. These 
parton-to-photon fragmentation functions therefore either have to be modeled in 
some way and/or constrained using experimental data. 

The fragmentation component represents a fraction of the inclusive prompt pho- 
ton signal which grows with the center-of-mass energy of the collision. While it re- 
mains subleading at fixed target energies, it becomes dominant at collider energies. 
On the other hand, most collider experiments - apart from the PHENIX experiment 
at RHIC [in], but, in particular, the TeV collider experiments CDF and DO at the 
Tevatron, ATLAS and CMS at the LHC - do not measure inclusive photons, be- 
cause at these energies the inclusive prompt photon signal would be swamped by a 
large background of secondary photons from decays of fast neutral mesons (mainly 
TT*^, as well as 77, etc.). Instead these experiments impose isolation criteria on the 
hadronic final states of photon candidate events, requiring that the photon be not 
accompanied by more than a prescribed amount of hadronic transverse energy in 
some given cone about the photon. An analogouil] criterion can be implemented in 
parton level calculations. The isolation cuts do not only suppress the background, 
they also substantially reduce the (F) component. Yet some fraction of the (F) com- 
ponent may survive and affect shapes of various tails of distributions, especially for 
correlation observables. Of course the sensitivity is even larger when loose isolation 
cuts are applied. In this article we wish to stress the interest of photon-jet correla- 

^The experimental request may also impose a veto on charged tracks in the vicinity of the pho- 
tons. However such vetoes cannot be implemented in partonic level calculations: a full description 
of the hadronised final state would be necessary. 
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tion observables, in particular in constraining the photon fragmentation functions, 
which requires to go beyond the lowest order. Similar studies have been performed 
in ref. [H]. However in these works the (F) component was calculated at Lowest 
Order (LO) only. 

In a previous article [12] we proposed a critical re-examination of the status 
of single prompt photon production in hadronic collisions in the light of recent 
experiments, which was based on a Next-to-Leading Order (NLO) calculation of 
both (D) and (F) provided in the form of a partonic Monte-Carlo code, JETPHOX 
|13j . The present paper aims at supplementing this previous work with a study 
of photon-jet correlations using the same tool, for the presentation of which we 
refer to [12]. The article is organised as follows. In section [21 we examine the 
magnitude of the fragmentation component on the photon-jet angular distribution 
at the Tevatron. In section [3] we discuss the jet rapidity distribution in photon + jet 
associated production, as a possible way to help constrain the uncertainties on the 
gluon distribution function. In sect. IHwe then discuss the potential of photon-jet 
correlations measured at RHIC without isolation as a tool to constrain the photon 
fragmentation functions. Finally, section [5] gathers our conclusions. A similar study 
dedicated to fragmentation into hadrons will be discussed in a future article. 



2 Photon-jet angular distribution 

An observable expected to receive a distinctive contribution from the (F) compo- 
nent is the photon-jet angular distribution which has been measured by the CDF 
Collaboration [TU [15] and is defined as follows. 

At lowest order (LO), corresponding to 2 2 kinematics, cos^* is the cosine 
of the angle between the photon direction and the beam axis in the center-of-mass 
system of the partonic subprocess. It also coincides with cos 6** = tanhy* where 
y* = {y^ — ?/jet)/2. This angular distribution is expected to receive a dominant 
contribution from the (F) component when cos^* becomes close to 1. Indeed, at 
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lowest order, the (D) component proceeds via a t-channel quark exchange yielding 
a behaviour ~ 1/(1 — cos^^*) for the partonic amplitude squared, whereas the (F) 
component involves also gluon exchange in the t-channel, yielding a behaviour ~ 
1/(1 — cos 6**)^. On this ground, one thus expects (F) to take over for cos 6'* values 
close enough to 1. 
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Figure 1: Sensitivity of the distribution of cos6'* in photon + jet at NLO to the 
isolation parameter e = E't max /Pt- Top left: Direct (D) contribution only. Top right: 
Fragmentation (F) contribution only. Bottom left: total (D) + (F) contribution. 
Bottom right: differences in (D), (F) and total (D) + (F) between e = 0.3 and 0.05. 

Depending on the cuts applied, the cos 9* dependences coming from the partonic 
transition matrix element squared may be blurred by an extra dependence com- 
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ing through the parton luminosity. Focusing on the direct (D) contribution, and 
parametrizing the LO phase space as 

diPr)'^ dy-r dy^et = ^ d{p*)^ dys d cos 9* (2.2) 

in terms of the variables cos6'* = tanhy*, yB = (1/7 + l/jet)/2 and p* = pj^ coshy*, 
the LO distribution in cos 6* reads: 

-J—a^ = 11 I ^Vb dp* Giip{xi) Gj/p{xj)- ^^^^ ^ (2.3) 

dcosu* J ' Ji J dcosO* dys dp* 

with 

x,, = ^ e^y- (2.4) 

In particular, one of the two parton distribution functions (pdfs) involved in the 
distribution in cos^* has an argument x which grows with y* i.e. with cos 6'* at 
fixed pj', so that this pdf decreases (towards zero) if cos^* increases (towards one). 
In absence of extra cuts, this decrease actually takes over the growth of the partonic 
cross section with growing cos 6* over the whole range. At LO, this can be neutralised 
by imposing cuts on ys and p* independent from y* [HI [15], so that the integration 
over ys and p* in eq. (12.31) yields cos 6** independent factors. Similar procedure 
and conclusion hold for the fragmentation (F) component, at least in the absence of 
isolation. 

Beyond LO the definition of cos 6'* has to be extended. This extension is not 
unique, and various definitions can be found in the literature. Here we take 

COS0* = tanhy* (2.5) 

where 

y* = \ (y-r - ?/leadingjet) , (2.6) 



^An alternative possibility is the one used in [H,, which, in short, combines the several jets of 
a multijet final state into one so-called superjet recoiling against the photon, in order to stick to a 
2^2 kinematics as close as possible. See [13] for more details. 
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y-y — yioadingjot being the differencci between the rapidity of the photon and the ra- 
pidity of the leading jet, i.e. the jet of highest transverse energy. Furthermore, 
the higher order contributions to the angular distribution involve an extra convo- 
lution smearing over kinematical configurations, so that the interpretation of this 
distribution beyond LO is less transparent. 

Besides, measurements at colliders most often involve isolated photons, in which 
case the (F) component is quite reduced. Namely, when the hadronic transverse 
energy accompanying the photon is required to be smaller than Exmax, the (F) 
contribution is roughly proportional to (1 — Zc) — ETmax/Pri the width of the 
support [ 2;^ , 1 ] of the convolution with the photon fragmentation functions. If 
Ermax is chosen such that this ratio is always small, the dominance of the t-channel 
gluon exchange from the (F) component is never effective; yet one might still expect 
a sizable distortion of the angular distribution for cos^* close enough to 1. We note 
also that if Exmax does not scale with pj., the isolation constraint implied on the 
fragmentation variable z, z > Zc with 

1 max 1 ^' 9 7^ 

~ p* (1 -COs2 0*)l/2 ^ • > 

induces, through the convolution over z^ an extra dependence on cos 9* which con- 
torts (amplifies somewhat) the growth of the (F) contribution provided by the par- 
tonic transition matrix elements alone. 

In a preliminary study [15] subsequent to the analysis published in [Ej and based 
on a data set with larger statistics and extended towards lower values of pj^, the CDF 
Collaboration found a discrepancy between the measured cos 9* distribution and the 
theoretical prediction of pjj. This subsequent preliminary CDF analysis concluded 
that extra dijet-like contributions involving t-channel gluon exchange would be nec- 
essary to bridge the gap, and that these extra contributions might come from NLO 
contributions to the (F) component. Since the prediction of [llj involves an account 
of (F) at LO only, we have revisited this observable and computed the effects of 



■^Note that the definition of cosO* given by (|2.5p . (|2.6p refers to a quantity which is invariant 
under longitudinal boosts along the beam axis. 
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accounting for (F) at NLO. The CDF Collaboration adopted a procedure to patch 
together the contributions from data corresponding to two regions that were distinct 
in p* and t/B though overlapping in cos^*, in order to maximise the range in cos^* 
displayed on one and the same plot. In particular the distribution was normalized 
to 1 in the bin farthest from cos^^* = 1, and the data sets from the two regions 
were normalized to each other in one overlapping bin in cos^^*. We understand this 
approximate procedure to have been dictated by the use of limited statistics and 
precision of Run I data but we did not follow the same procedure in our study for 
several reasons. Firstly, the normalization to 1 in the bin farthest from cos 6'* = 1 
aims at getting rid of the numerical factor coming from the integration of partonic 
luminosity. This is fine as long as only one partonic subprocess contributes - or at 
least, when one yields a much greater contribution than all the others. However, in 
the present case, the distribution is of the form 

E^'"^ (2.8) 



d cos 9* g d cos 9* 

i.e. a linear combination of contributions coming from several subprocesses s. In 
particular, considering (D) only, gq (or gq) initiated and qq initiated processes which 
contribute at LO have distinct functional dependences on cos 6**. Secondly, the 
integrated partonic luminosity factors depend not only on the subprocess s 
but also on the integration regions in phase space. The relative weights of the 
subprocesses from the contributions (D) and (F) thus differ in the two regions defined 
by CDF as "1" and "2" . The normalization enforced by the matching procedure of 
CDF is not harmless on the impact of the NLO correction to the (F) contribution, 
and might reduce the impact of this correction. 

Therefore we did not stick to the CDF study. We focused on a study of the 
magnitude of the (F) contribution, without making a direct comparison with the 
CDF Run I data. Our study has been made for =1.96 TeV, with the following 
definitions and kinematic cuts: \y^\ < 0.9, px > 30 GeV, p^^°-'^^^9 ^"^^ > 25 GeV. The 
jets were defined according to the DO midpoint algorithm [22], with cone aperture 
Rc = 0.7. The photon isolation required that in a cone of aperture R = 0.4 in 
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rapidity and azimutal angle around the photon direction, the fraction of maximal 
hadronic transverse energy ETmax/Pr be less than a prescribed value e, which we 
varied from 0.05 to 0.3. The further cut 45 GeV < p* < 55 GeV was imposed. 
The pdf set CTEQ 6.1 was used together with the BFG set II for the fragmentation 
functions, with the scale choice /i = Mp = M = pj,/2. 

We have considered three ingredients which may affect the size of the contribution 
of the (F) component. One is the account for the NLO corrections to the many 
subprocesses; another one concerns the uncertainty on the fragmentation functions; 
yet another one deals with a possible mismatch between the implementation of 
isolation at the partonic vs. hadronic level. 

Let us first consider the impact of the NLO corrections to the (F) component. 
For the standard scale choice /i = Mp = M = the effect is to multiply the 

component (F) at LO by about a factor two. From the top-right fig. [H the impact 
on the total (D) + (F) both at NLO amounts to an increase by 4 % in the upper 
cos 6** range w.r.t. (D) at NLO + (F) only at LO for e = 0.05. 

Is it possible to increase the (F) contribution by modifying the fragmentation 
functions? When a stringent isolation cut is required on the photon candidates as 
in the CDF experiment, the (F) contribution involves the photon fragmentation 
function at 2; > i-e. rather close to 1. In this region the fragmentation functions 
are dominated by their so-called anomalous parts predicted by perturbative QCD. 
Their poorly known non perturbative parts, which would be the only adjustable 
ingredients, play no role: thus the (F) contribution to the cos^* distribution is 
rather tightly constrained. 

We have tackled the issue of the account of isolation at the partonic vs. hadronic 
level by varying the value of the isolation parameter e from 0.05 to 0.3. As already 
mentioned in table 2 of [16] in the case of the inclusive cross section and as can be seen 
on fig. [U the separate (D) and (F) contributions do depend strongly on e at NLO; 
yet strong cancellations turn out to occur between (D) and (F) so that the total (D) 
-|- (F) depends on e only very mildly, at least as long as the infrared sensitive term 
as In e in (D) does not become large - otherwise the fixed order prediction becomes 
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unreliable. Therefore, changing the calorimetric isolation parameter by as much 
as a factor of six does not modify the total contribution to the cos^* distribution 
significantly. 

To summarize, the results of our calculations show that the idea of playing with 
the fragmentation component as suggested in the CDF analysis turns out to be 
ineffective in the conditions which we have considered. It would be worthwhile to 
perform a quantitative analysis of the much larger statistics data set gathered in 
Run II, without relying on the questionable matching procedure used in the CDF 
Run I analysis. 

Let us recall that, beyond the isolation requirements and the refined analysis to 
improve background rejection, the CDF measurement of isolated photons required 
the statistical subtraction of a contamination of photon candidates coming from 
neutral hadrons - mainly 7r°, plus t], p^, u etc. The background was removed statis- 
tically by exploitation of the expected difference in the resulting shower profiles or by 
the different conversion probabilities using Monte Carlo simulations [17] . The latter 
relied on hadronisation models suited to describe the bulk of hadronisation. On 
the other hand, the very small fraction of hadronic events which pass the isolation 
cuts corresponds to the tail of fragmentation at large z which is not constrained by 
the data. These background events yield namely dijet-type contributions involving 
t-channel gluon exchange, which might explain part of the discrepancy observed, 
and the distribution in cos 6'* at cos 6'* 1 might provide an enhanced sensitivity 
to this contamination w.r.t. other prompt photon observables. 

Let us add a comment on the comparison with the situation for the photon- 
photon azimuthal angle distribution in photon pair production [18]. In the latter 
case, some of the higher order contributions involving one direct photon and one 
photon from fragmentation provided a coUinear logarithmic enhancement at low 
azimuthal angle, when a hard jet recoils against the photon pair at low relative 
angle. No such phenomenon occurs in the photon-jet case, since the jet considered 
in that observable is always the leading jet of the event which roughly recoils against 
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the photon. 



3 Jet rapidity distribution 

Despite years of intense work, a proper understanding of the uncertainties on the 
gluon distribution function is still lacking [19j: due to shape assumptions, one can 
find some regions with small errors despite the lack of data points. More precisely, 
as shown by [20], uncertainties on the gluon distribution at very low x {x < 10~^, 
where there are no constraining experimental data points) obtained by CTEQ6.5 
and MRST2001e do not overlap. More flexible shapes modeled with neural networks 
and fitted using a genetic algorithm from NNPDF [21] give much larger error bands 
overlapping with the ones from CTEQ and MRST. We note that dynamical PDFs 
generated radiatively from valence like input at low scales may be another approach 
which yields smaller uncertainties, [22], see also fig. 3 of [23]. 




Figure 2: Distribution of leading jet rapidity in photon + jet associated production, 
at y-y = for various pdf sets and scale choices. 
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Figure 3: Distribution of leading jet rapidity in photon + jet associated production, 
at = 2.5 for various pdf sets and scale choices. 

Is it possible to be less sensitive to shape assumptions by using photon+jet cor- 
relations? The Tevatron experiments CDF and DO both during Runs I and II, have 
been measuring photon-jet correlations, which explore the short distance dynamics 
in a more constrained way than inclusive photon production. A recent comparison 
between DO data and JETPHOX has been performed j24j for the distribution da/dy^ 
vs. pj.. An interesting studjfl [25] of the possibilities of the CMS experiment on 
this distribution has recently appeared. Among the other correlations which can 
be studied, let us mention the distribution of jet rapidity at fixed photon rapidity, 
integrated over the photon transverse momenta above some Prmin- large rapidi- 
ties, the main contribution comes from the subprocess qg g7 (or qg 57) where 
the initial state gluon is at quite low x, down to (9(a few 10~^) while the x of the 
initial state (anti)quark is ~ (9(10^^). This correlation observable is thus sensitive 
to gluons in a low x region overlapping with the one explored at HERA. 
^Yet we notice that this study accounts for the fragmentation contribution at LO only. 



12 



Figs. [2] and [3] show the JETPHOX predictions for various pdf sets for the distri- 
bution of jet rapidity, for the photon rapidities = and 2.5, respectively at the 
Tevatron for = 1.96 TeV. The jets are defined according to the DO midpoint 
algorithm with cone aperture Rc = 0.7 [26], and the discussion at NLO refers to the 
leading jet i.e. the jet with highest px- The cross section is integrated over photon 
transverse momenta larger than 30 GeV and over jet pt larger than 20 GeV. The 
choice of scales is = M = Mp = p'r/'^- Besides the prediction usingj the pdf 
sets CTEQ 6.6 [27] and MSTW08 [28] resulting from global fits, we also show the 
prediction with the BBS set [29], an example of a set modeled through dynamical 
generation, which has a quite different gluon pdf in the low x region, to illustrate the 
sensitivity of this observable to low x gluons. However we cannot draw any definite 
conclusions from this observable alone since the dependence of these predictions on 
the scale choice at NLO is as large as the spread with respect to the various pdfs 
used. An error analysis taking into account the detailed information provided by 
MSTW08 is beyond the scope of this paper. 

4 The Photon Fragmentation Function 



The large-pT photon-jet correlations also give access to the Photon Fragmen- 
tation Function (PFF) which has rarely been measured. Actually only two LEP 
experiments, ALEPH [31j and OPAL [32] measured the PFF. However it is difficult 
to observe a photon in a large hadronic background and hence the PFF has mainly 
been measured for large values oi z = lE^j ^/s. Good agreement is found between 
these data and two NLO theoretical results [30| 133]. 

The hadroproduction of large-p^ photons and jets should also allow to measure 

the PFF. First let us consider the direct subprocess qg q'y in which the final 

^The photon fragmentation functions (PFF) are from BFG set II [30]. Due to the DO isolation 
requirement - less than 2 GeV of accompanying hadronic transverse energy in a cone of radius 0.4 
in azimut and rapidity around the photon direction - the prediction depends only marginally on 
the PFF choice. 
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Figure 4: Cross section da/dz-y for the BFG set II. 

photon transverse momentum p7p is balanced by the g-jet transverse momentum 
such that = — t^j^tj = 1- The PFF is not involved in the description of 

WPt II 

this reaction. The PFF manifests itself, for instance, in the subprocess reaction 
gq —>■ gq followed by the gluon ((? — * 7 + X) or the quark {q —>■ 'J + X) coUinear 
fragmentations described by the distributions D^/g{z^, Mp) and D^/q{z^, Mp), the 
large scale Mp being of order (pt)"^- At Leading Order (LO) we have z^ < 1 and 
the cross section da^^^^/dz^ is directly proportional to the functions D^/a{z^, Mp) 
(a = q,g). Therefore at leading order the total cross section da/dz^ is given by the 
sum of the direct contribution (proportional to 5(1 — z^)) and of the fragmentation 
contribution containing the PFF D^/aiz-y, Mp). Contrarily to e^e~ experiments we 
have to include a direct contribution in the cross section and to stay away from 
z^ = 1 to increase the sensitivity to the PFF. 

When HO corrections are taken into account, more jets can be present in the 
final state and z-y may be larger than one. With three partons in the final state, z^ 
can be different from one also in the direct contribution. 

The variable depends on the two transverse momenta and pp^, none of them 
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Figure 5: da/dz^ for tlie two BFG sets I and II. 

being a priori fixed for a given value of z^. However is directly related (at least at 
LO) to the parton momentum involved in the hard subprocess. Therefore, if we vary 
jTy* to obtain different values of z^ (keeping p^. fixed), the theoretical cross section 
da /dz^ will reflect the z.y dependence of D^/a{z^, Mp) and the p'^^ dependence of the 
subprocess, thus blurring the 2;^-dependence of the PFF. On the contrary, if we keep 
j?^;* fixed and vary pj,, we obtain a z^-dependence of da/dzj coming dominantly 
from the fragmentation function. Therefore we propose to measure the PFF in 
experiments in which the jet momentum is kept fixed and the photon momentum is 
varied. 

Let us finally note that the observed photon must not be isolated, which would 
considerably reduce the fragmentation contribution. This possibility exists when the 
photon-pT is not too large, as it is the case at RHIC for pj, < 16 GeV [TO t IM j 135 1 136] . 
Therefore we choose the photon and jet momenta in agreement with the experiment 
performed at RHIC, i.e. in our numerical analysis we use 3 < pj. < 16 GeV, 11 GeV 
< < 13 GeV and ^/s = 200 GeV. For the rapidities we take -| < < | 
and — 1 < yj'^* < 1. We use the CTEQ6M parton distribution functions [TTj. and 
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C = 1/2 
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C = 2 



0.2 0.4 0.6 



Figure 6: Scale dependence of da/dz^ for a common variation of scale fj, = M = Mp. 
renormalization and factorization scales sriven by C ■ with \ <C <2. 

Discussion of numerical results 

In Fig. m we present the results for the cross section da/dz^, calculated with the 
scales /i = M = Mp = p'p and the set 11 of the BFG fragmentation functions [30] . 

The jet is defined by the midpoint algorithm [26] with -Rcone = 0.7. We ob- 
serve that at small values of z^, the fragmentation contribution is much larger than 
the direct one. Although at next-to-leading order, the fragmentation and direct 
contributions cannot be considered as independent physical channels, this feature 
persists for other scale choices. Therefore the photon fragmentation functions should 
be measurable in RHIC experiments at small values of 2;^, a region which was not 
accessible to LEP experiments. 

In Fig. [5] we compare the predictions obtained with the sets 1 and 11 of the 
BFG fragmentation functions. The gluon fragmentation function of set 1 is smaller 
than the one of set 11 at small z^, leading to a significant difference in the da/dz^ 
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Figure 7: Separate renormalization scale (yu) dependence of da/dz^ for fixed common 
factorization scales M = Mp = p'^ ■ 



distribution, as can be seen from Fig. [51 This result shows the sensitivity of this 
reaction to the photon fragmentation functions. 

We find that the scale dependence at NLO is about ±20% as can be seen from 
Fig. [Hi There are three scales involved: the renormalisation scale /i, the initial state 
factorisation scale M (we set the factorisation scales for both initial state particles 
equal to M), and the final state factorisation scale Mp. The dependence of the cross 
section on these scales is quite different. 

In Fig. [6] we vary all three scales simultaneously. If we vary the final state 
factorisation scale Mp only, keeping /i = M fixed, the fragmentation component 
increases logarithmically with Mp^ while the direct contribution decreases, such that 
the leading logarithmic dependence cancels in the sum and the overall dependence 
on Mp is rather weak. We also checked that the dependence on the initial state 
factorisation scale is very weak, such that the overall scale uncertainty is dominated 
by the renormalisation scale dependence. This is demonstrated in Fig. [3, where 
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BFG set II, M=Mr=pT'V2 



M=Pt'V3 
M=Pt'*'/2 



Figure 8: Region of larger stability w.r.t. renormalization scale dependence. 



we vary the renormalisation scale fj, only, keeping M = Mp fixed to pj^. However, 
we can find a region where the cross section is more stable against variations of fi, 
which is in the vicinity of M = Mp = and p't'/4: < jj < Pj;*/2, as can be seen 

from Fig. [HI This optimal behavior is obtained in the small- 2; domain in which the 
fragmentation contribution is large. 

Fig. [9] shows that the size of the higher order corrections is much larger for the 
fragmentation component than for the direct component. 

5 Conclusions 

In this article we have studied photon-jet correlations in hadronic collisions based on 
the NLO program JETPHOX, which is a Monte Carlo program of partonic event gen- 
erator type which incorporates NLO corrections to both direct photons and photons 
from fragmentation. The program is flexible to account for user-defined kinematic 
cuts and photon isolation parameters. It is available at the following web site [13] . 
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frag LO 



direct NLO 
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Figure 9: Respective sizes of HO corrections in fragmentation and direct components 
for a given common scale choice. 

Correlation observables offer in general a larger sensitivity to the short distance dy- 
namics than one-particle inclusive observables. We studied the photon-jet angular 
distribution cos^* and the jet rapidity distribution in view of possible constraints 
on the parton distribution functions in the proton, in particular the gluon. 

Furthermore, correlations involving unisolated photons such as the ones mea- 
sured at RHIC provide a means to constrain the photon fragmentation function in 
a region which was barely accessible by the LEP experiments. We study the ob- 
servable da/dz^, where z-^ can be reconstructed from the photon and jet transverse 
momenta. We argue that the z^ dependence of da/dz^ is coming dominantly from 
the fragmentation if we keep p^;* fixed and vary pj., and therefore propose to mea- 
sure the fragmentation functions at fixed pj^*. The NLO predictions still have a non 
negligible dependence on the renormalisation scale choice, which is related to the 
fact that the higher order corrections to the fragmentation component are large, 
exceeding 40%. 
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The program JETPHOX can also predict the tail of the distribution of transverse 
momentum Qt of photon+jet pairs at NLO accuracy for large enough Qt for the 
LHC, as well as any other correlation insensitive to multiple soft gluon emission, it 
thus can also help to normalize the Monte Carlo event generators suited to describe 
the lower Qt range of the distribution of transverse momentum of photon+jet pairs 
and other less inclusive observables. 

Note added 

While we were completing the present article, we became aware of a work by Stavreva 
and Owens [38j , similar to the one reported here, but focusing on charm and bottom 
jets in the final state. 
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